Generation of coherent acoustic phonons in piezoelectric semiconductor 
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We review some experimental and theoretical aspects of coherent acoustic phonon generation in 
piezoelectric semiconductor multiple quantum wells. In order to model more advanced and compli- 
cated nano-acoustic devices, a macroscopic continuum theory for the generation and propagation of 
coherent acoustic phonons in piezoelectric semiconductor heterostructures is presented. The macro- 
scopic approach is applicable in the coherent regime, and can be easily utilized to study coherent 
acoustic devices based on piezoelectric semiconductor heterosutructures. For each phonon mode, 
the corresponding coherent acoustic field obeys a loaded string equation. The driven force has 
contributions from the piezoelectric and deformation potential couplings. We applied the theory 
to model the generation of coherent longitudinal acoustic phonons in (OOOl)-oriented InGaN/GaN 
multiple quantum wells. The numerical results are in good agreement with the experimental ones. 
By using the macroscopic theory, we also investigated the crystal-orientation effects on the genera- 
tion of coherent acoustic phonons in wurtzite multiple quantum wells. It was found that coherent 
transverse acoustic phonons dominate the generation for certain orientation angles. 

PACS numbers: 62.25.+g, 63.22.+m, 43.35.+d, 78.47,+p 
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I. INTRODUCTION 

Ultrafast phenomena in condensed matters has been 
for many years the subject of intense experimental and 
theoretical activities. Progress in femtosecond lasers and 
ultrafast spectroscopy has enabled us to investigate the 
initial relaxation of the nonequilibrium photoexcited sys- 
tems Typical examples of time-resolved techniques 
include pump-probe measurement 0, |j| , transient four- 
wave-mixing p| , and time-resolved fluorescence Q ■ 

Optical excitation in semiconductors creates nonequi- 
librium electron and hole distributions as well as coherent 
interband and intraband polarizations. Time evolution 
of these coherent electronic polarizations has led to the 
observation of different coherent phenomena and oscilla- 
tions including Bloch oscillations in superlattices 0, Fj , 
heavy and light hole quantum beats in quantum wells |8| , 
Rabi flopping in semiconductors |9j , coherent dynamics of 
excitonic wave packets , wave packet oscillations , 
THz emission from asymmetric double quantum wells 
|12| , and far-infrared emission from asymmetric quantum 
wells . At a time scale of a hundred femtoseconds to 
a few picoseconds, another macroscopic coherence phe- 
nomenon, the coherent phonon oscillations, is observed. 
While the macroscopic coherence quantities correspond- 
ing to the electronic system are intraband and interband 
polarizations which contribute to the above mentioned 
electronic oscillation phenomena, coherent phonon oscil- 
lations result from the existence of quantum averages 
of the phonon creation and annihilation operators 



Time-resolved observations of coherent optical phonons 

have been reported for semiconductors, e.g. GaAs [l5| 

and Ge 0, semimetals, e.g. Bi and Sb |17|. cuprate 
superconductors |l8j . and a number of other materials 
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On the other hand, generation and detection of co- 
herent acoustic phonons has also been demonstrated in 
bulk materials [20L l2l| , as well as artificial heterostruc- 
tures, e.g. superlattices H3,|2^|, quantum dots [24J, and 
nano-particles (2{|. Recently, Sun et al. have observed 
large amplitude longitudinal-acoustic (LA) phonon os- 
cillations in piezoelectric InGaN/GaN multiple quantum 
wells (MQW) 26]. In addition to the electron-phonon 
deformation potential coupling mechanism, the genera- 
tion of coherent LA phonons within the MQW structure 
was found to be dominated by a strong piezoelectric cou- 
pling mechanism. This is due to the large GaN piezo- 
electric coefficients and the strong built-in piezoelectric 
fields in the strained InGaN epilayers. Coherent con- 
trol of the LA phonons within a few oscillation cycles 
was also demonstrated in this system [27|, H^. Such a 
nitride-based MQW structure acting as a coherent THz 
phonon source opens a new passage to phonon engineer- 
ing, which might significantly enhance the performance 
of nano-scale solid-state devices. 

In this paper, we review some experimental and the- 
oretical aspects of coherent LA phonons in piezoelectric 
MQWs. We also present a macroscopic theory based on 
continuum elastic dynamics and electrodynamics. This 
macroscopic approach can be easily extended to analyze 
more complicated coherent acoustic nano-devices. We 
shall first review some experimental results in section 2. 
The microscopic theory of coherent LA phonons in In- 
GaN/GaN MQWs is discussed in section 3. In section 4, 
we present our macroscopic theory. We solved the loaded 
string equation and introduce the sensitivity function in 
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Section 5. In Section 6 we apply the macroscopic ap- 
proach to investigate the crystal orientation effects on 
coherent acoustic phonon generations. Finally we make 
a conclusion in Section 7. 



II. TIME-RESOLVED PUMP-PROBE 
MEASUREMENT OF COHERENT LA PHONONS 

Femtosecond lasers have proven to be powerful tools 
for studying the dynamical behavior of photoexcited elec- 
trons and holes in semiconductors with a time resolution 
as short as a few femtoseconds. Here we review some 
results of our pump-probe measurements on coherent 
acoustic phonons in GaN heterostructures. FigQ]shows a 
schematic diagram of the transient transmission measure- 
ment setup, which is usually referred to as a pump-probe 
setup. The output of the femtosecond laser is divided into 
two beams; one is used as the pump beam and the other, 
much weaker, beam is used as the probe. A mechan- 
ical delay stage introduces a time delay between these 
two synchronized pulse trains. The pump pulse photoex- 
cites the semiconductor and modifies its optical proper- 
ties, e.g., absorption. The excitation relaxes within a few 
tens of femtoseconds to several hundred picoseconds. By 
measuring the induced transmission change of the probe 
pulse as a function of time delay, we can study the time- 
resolved carrier dynamics of the photoexcited semicon- 
ductors. 

We first discuss the basic mechanisms responsible for 
photoexcitation and detection of large amplitude coher- 
ent LA phonons in strained InGaN/GaN MQWs with 
built-in piezoelectric fields. The UV femtosecond pulse 
first photoexcited electrons and holes within the MQWs. 
These photogenerated carriers are spatially separated 
due to the strong built-in piezoelectric field. The photo- 
carriers in turn screen the piezoelectric field and change 
the equilibrium state of the lattice through piezoelec- 
tric coupling. The lattice relaxes toward the new in- 
homogeneous configuration, and a displacive coherent 
phonon oscillation is thus initiated. The induced acous- 
tic phonon oscillations result in piezoelectric field modu- 
lations, which cause variations in the absorption through 
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FIG. 1: Experimental setup of UV femtosecond pump-probe 
measurement. 



the quantum confined Franz-Keldysh (QCFK) effect p^ . 
After photoexcitation, the observed probe transmission 
changes decay with time as the propagating acoustic 
phonons leave the MQW, resulting in a decay time con- 
stant proportional to the total number of wells. It is 
important to note that the decay of the transmission is 
related to the phonons leaving the MQW and not the 
actual decay of the phonon modes. 

FigUl shows the measured probe transmission changes 
as a function of probe delay for a 14-period 50 A/43 A In- 
GaN/GaN MQW sample at wavelengths of 390 nm (3.177 
eV) and 365 nm (3.39 eV). The average incident pump 
power was 20 mW. And the average 2D/3D photocar- 
rier densities were 9 x 10 12 cm~ 2 /1.8 x 10 19 cm -3 and 
1.5 x 10 13 cm- 2 /3 x 10 19 cm~ 3 for the 390 and 365 nm 
traces, respectively. After the pump excited carriers and 
caused a large transient transmission increases at zero 
time delay, a cosine-like transmission oscillation could be 
observed on top of the carrier cooling background signal. 
The amplitude of the cosine-like transmission modula- 
tion, AT/T , was on the order of 10~ 2 . Within our ex- 
perimental resolution, the observed oscillation frequency 
was found to be independent of the pump/probe pho- 
toenergy or pump fluence. The transmission oscillations 
can be fit using cosine functions with phases of zero or 
7r, which suggests that the oscillations are displacive in 
nature [l7j. The observed cosine- like oscillation is con- 
sistent with the picture that a new equilibrium configu- 
ration for the lattice system is set up by photo carrier 
screening of the strain-induced piezoelectric field. We 
shall discuss the generation mechanism and the displacive 
property of coherent acoustic phonons in more detail in 
the following sections. 



III. MICROSCOPIC THEORY OF COHERENT 
ACOUSTIC PHONONS 

Recently, Sanders et al. presented a microscopic the- 
ory for the generation and propagation of coherent LA 
phonons in strained InGaN/GaN MQWs [U- The 
authors used a density matrix formalism to treat the 
photoexcited electronic system and coherent acoustic 
phonons. And it was found that under typical experi- 
mental conditions, the microscopic theory can be simpli- 
fied and mapped onto a loaded-string problem that can 
be easily solved. Here we review some important results 
of the microscopic theory for comparison with the macro- 
scopic approach. 

The heterostructure is assumed to have a cylindrical 
symmetry along the crystal c-axis which is along the z 
axis, i.e. the (0001) direction. Thus the quantized state 
of electrons has a localized wavefunction along the z axis 
due to quantum confinement of the heterostructure, while 
in the x — y plane, the electrons are free to move. For 
simplicity, we also neglect the discontinuities of the elas- 
tic constants within the heterostructure as a first order 
approximation. The phonon eigenstates are just bulk- 
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FIG. 2: Measured transient transmission changes versus 
probe delay for a 14-period 50A/43A InGaN/GaN MQW 
sample with an average pump power of 20 mW. (a) Laser 
wavelength was 390 nm. (b) Laser wavelength was 365 nm. 
Cosine-like oscillation can be observed on top of the back- 
ground signals. Please notice the n phase shift between traces 
(a) and (b). Oscillation dephasing time was on the order of 8 
ps. 



indices, k is the in-plane wave vector for electrons in the 
subbands. The interaction matrix elements describing 
deformation potential and piezoelectric scattering are 



M a , = 



h 2 



2p (Hu; q )V 

|e|e 3 3 



£oo£ s (<7) 



iqV« in ,(k,q)- 



(3) 



where V is the crystal volume and po is the density. The 
first term in Eq. describes deformation potential scat- 
tering while the second term describes screened piezo- 
electric scattering. The factors V and V are related to 
the electron confinement wavefunctions along the z axis. 

In order to apply the generalized density matrix for- 
malism, one define the electron density matrix 



(4) 



where ( ) denotes the quantum and statistical average 
of the nonequilibrium state of the system. The inter- 
band components of the density matrix, N°', (k, t) and 
7V^; c n (k, t), describe the coherence between conduction 
and valence electrons in subbands n and n' and are re- 
lated to the optical polarization. The intraband com- 
ponents of the density matrix N"'", (k, t) describe cor- 
relations between different subbands of the same carrier 
type if n ^ n'. Hn = n', N%£(k, t) = f%(k, t) is just the 
carrier distribution function for electrons in the subband 
state. 

The coherent phonon amplitude of mode q is defined 
to be El 



D q (t) = (bl(t)+b- q (t)}- 



(5) 



The coherent phonon amplitude is related to the macro- 
scopic lattice displacement u(z, t) through the relations 



like plane-waves. Due to the cylindrical symmetry, only 
q = qz acoustic phonons are coupled by the electron- 
phonon interaction. The free LA phonon Hamiltonian 
can be written as 



(1) 



where b q and b q are creation and annihilation operators 
for LA phonons with wave vector q = qz. The phonon 
dispersion relation is simply u> q = c s \q\, where c s is the 
LA phonon sound velocity for propagation parallel to z. 

The interaction Hamiltonian for LA phonon and elec- 
trons is 

7~LeA — M% >n ,(k,q)(b q + bl)cl n ^c a , n , tk , (2) 

ct,n,n' ,k,g 

where c and c* are the creation and annihilation oper- 
ators for electrons. The index a — {c, v} refers to con- 
duction or valence subbands, and n, n' are the subband 



2 Po (h^ q )V 



(6) 



The equation of motion for these statistical operators 
can be obtained using the Heisenberg equation: dA/dt = 
((i/h)[H, A]). The electrons also couple to the laser field 
through dipole interaction. The coherent phonon ampli- 
tudes satisfy the driven harmonic oscillator equations 



d 2 D q {t) 
dt 2 



u 2 q D q {t) 



2LO q s~\ 

a.n.n' .k 



xTO(M) (7) 



subject to the initial conditions: D q (t = — oo) = dD q (t = 
—oo)/dt = 0. The equations of motion for coherent LA 
phonons are coupled to the electronic intraband polar- 
ization TV"'", (k, t) whose dynamical equations can also 
be derived from the Heisenberg equation. Based on a 
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linear dispersion relation for LA phonons, one can trans- 
form the above dynamic equation for D q (t) into a loaded 
string equation for u(z,t) 

d 2 u{z,t) j d 2 u(z,t) _ 

~W ° s dz* ~ b( - z ^ W 

with the following source function S(z,t) 



a,n,n' k,g 



'2hc s \q\ 
PoV 



M a n ,Jk, q y 



(9) 



The initial conditions for u{z, t) are: u(z, t — — oo) = 
<9u(z, t = —oo)/dt — 0. In the following section, we shall 
derive the same loaded string equation from a macro- 
scopic approach with a simplified source function. 



IV. MACROSCOPIC CONTINUUM THEORY 
FOR GENERAL CRYSTAL ORIENTATIONS 

In this section, we develop a continuum elastic model 
for the generation of coherent acoustic waves in nitride- 
based nanostructures. The model is based on the macro- 
scopic constituitive equations taking into account both 
the piezoelectric and deformation potential couplings. 
The governing dynamical equations are the elastic wave 
equations coupled to Poisson equation. This approach is 
valid in the coherent regime since thermal and quantum 
fluctuations can be neglected. More specifically, these 
equations can be regarded as the corresponding opera- 
tor equations with appropriate quantum averages of the 
nonlinear terms. Here we consider coherent phonon gen- 
eration in arbitrarily orientated MQWs. In order to con- 
sider the crystallographic effect, we use (x,y,z) to de- 
note the primary crystallographic axes of wurtzite ni- 
trides. The z direction is along the crystal c-axis and 
n is along the crystal growth direction. From the sym- 
metry of wurtzite crystals, the macroscopic properties 
depend solely on the angle 8 between n and the c-axis. 
Thus we may let n||[M)zZ]. From the Miller-Bravais nota- 
tion [3l]], (hOil) — (hOhl), so that the polar angle 8 can 
be expressed as a function of indices h and I only with 
cos 8 = ul/ y/ 'Ah 2 /3 + u 2 l 2 . Here u = \fafc is the inter- 
nal structure parameter, with a and c the usual hexag- 
onal lattice parameters. We use £ = n ■ r to denote the 
coordinate along the growth direction. 

The free energy density of a piezoelectric semiconduc- 
tor is mm 



F(T,Ei,€jk) — Fq(T) + -CijkitijCki — —EijEiEj 

-eijkE-jejk + ^2 11 "■<>■,■' (10) 

v=e,h 

where we have taken the temperature T, electric field 
Ei, and strain ej k as independent thermodynamic vari- 
ables. In the above expression, Cyjy, e^, and e^-fc are 



the isothermal elastic stiffness, dielectric, and piezoelec- 
tric tensors, respectively. The index v runs over carrier 
species (electrons and holes), p v is the number density of 
carrier species v, and d V i is the corresponding deforma- 
tion potential. Please note the convention that repeated 
indices imply summation is used throughout the text. 
The last two terms in Ea. (|10|l . represent piezoelectric and 
deformation potential couplings, respectively. From the 
free energy density, we can derive the stress tensor and 
electric displacement according to cr^ = dF/dtij and 
Di = —dF/dEi, respectively. We obtain the following 
constitution equations for piezoelectric semiconductors 



o'ij = Ciju^u — ekijEk + 



^ ^ 3ij dyj Pi/. 



D, 



ZijEj 



■ eikl e ki- 



(11) 



(12) 



As discussed in the previous section, we neglect the elas- 
tic discontinuities here and regard the semiconductor het- 
erostructure as a continuous medium. For InGaN/GaN 
heterostructures, this approximation is plausible since 
the elastic properites of InN is close to those of GaN. 
The macroscopic equation of motion is 



Po- 



d 2 Ui d<Ji 



dxj 



(13) 



where Ui is the lattice displacement relative to the static 
equilibrium state in the absence of photoexcitation. 

The electric displacement D satisfies the Poisson equa- 
tion 

V-D(£,t) =p sc &t) = \e\[p h (t,t) -pefot)], (14) 

where p sc is the space charge density. From the cylin- 
drical symmetry of the carrier distribution which varies 
along the growth direction of the heterostructure, these 
dynamical variables are functions of the spatial variable 
£ only. In the quasi-static approximation, we assume the 
electric field is irrotational, E = E sc h. By substitut- 
ing the constitution equation into dynamic and Poisson 
equations, we obtain 



dE sc 
: 9£ 



(15) 



Po 



= r ?; 



dE„. 



in- - " <>; + (16) 



Here we have introduced the following effective constants 
along direction n: 



= CikljUkTil 

e-j = e k ijn k ni 



(17) 
(18) 
(19) 
(20) 
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Tij is the Christoffel tensor, e is the effective dielectric 
constant, ej is the effective piezoelectric tensor, and d„i is 
the effective deformation potential coefficient for species 
v. Note that no summation over index i is assumed in 
Eq.(20). By eliminating the electric field, we obtain the 
following loaded wave equation 



on this direction vanishes. The other two eigenmodes 
are in the plane spaned by c and n. Their polarizations 
are quasi-transverse and quasi-longitudinal. For orienta- 
tion along [0001], these two modes are exact longitudinal 
acoustic (LA) and transverse acoustic (TA) phonons. 



pa 



d 2 u t 
dt 2 



r, 



where Tij = Tij 



etej/e is the effective Christoffel tensor. 
The two driving forces are 



/t,piezo(C,*) = =Psc(£,t) = \e\4[pe(£,,t) ~ Ph(£,,t)] 

(22) 



/i,def(£,*) 



i/—e,h 



/i )P ie Z o + /i,def, (21) V. LOADED STRING EQUATION AND 

DETECTION OF COHERENT ACOUSTIC 
PHONONS 

The loaded string equation can be solved by Green's 
function method as discussed in Ref 30. Here we solved 
the loaded string equation for the case of a InGaN/GaN 
MQW with 4 wells. The In composition is 10% and the 
well and barrier widths are L w = 22 A and L b = 90 
A, respectively. In Fig|3|we show the time evolution (from 
(23) the bottom trace to the top trace) of the coherent phonon 
strain field. The strain field is s(z,t) = du(z,t)/dz The 



1 <9<e 



/i.piezo is the driving force due to piezoelectric coupling 
and /i,def is caused by deformation potential. From 
Ea. (|22l) . the piezoelectric force is proportional to the 
space charge density and is parallel to the effective piezo- 
electric constant ej. On the other hand, the deformation 
potential force is proportional to the spatial derivative of 
the carrier density and is parallel to the effective defor- 
mation potential coefficient d V i for each carrier species 
v. Since the deformation potential constants d V i are ap- 
proximately the same for different directions i = x,y,z, 
it can be seen from Eq.(20) and l123[l that the deforma- 
tion potential force approximately points to the crystal 
growth direction n. 

To solve the loaded wave equation l|21|) . we first find 
the eigenmodes of the Christoffel tensor. The coherent 
acoustic waves are then expressed as superpositions of 
these eigenmodes. The driving forces for each eigenmodc 
are then obtained by projecting the general forces onto 
the directions of the eigenmodes. The eigenmodes are 
obtained by solving the Christoffel equation 



r, 



W j 



c 2 pom 



(24) 



where c is the phase velocity of the eigenmode and Wi 
is the corresponding polarization of the lattice displace- 
ment. The Christoffel equation is an eigenvalue equa- 
tion with poC\ being the eigenvalue and w\i the eigen- 
polarization. A = 1,2,3 correspond to the three eigen- 
modes. The general solution of the acoustic field is 
Ui = Yl\=\ 2 3 u \w\i- And the mode amplitudes u\ sat- 
isfy the following loaded string equation 
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d 2 u x 
dt 2 



2 d 2 u x _ 1 

— c-i —rr^r — — \J A.piczo ™r J \,dcl) 
PO 



(25) 



One of the cigcn-polarizations is always perpendicular to 
the plane spanned by c-axis and the growth direction N. 
This mode is an exact transverse mode. However, it will 
not be excited since the projection of the driving forces 



FIG. 3: schematic diagram of the time evolution of the co- 
herent phonon strain field and the spatial distribution of sen- 
sitivity function. The strain field was derived from the loaded 
string equation. Also shown is the calculated oscillating traces 
using the sensitivity function and strain field solution on the 
left figure. 
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results shown in Fig[3|correspond to strain distribution at 
time delays t = mT/2, where m = 0, 1, 2, . . . and T is the 
fundamental oscillation period. After the acoustic wave 
trains leave the MQW region, a static carrier-induced 
strain exists in the MQW. 

The transmission changes of the probe pulse due to 
the existence of coherent acoustic pho nons can be de- 
scribed using a sensitivity function |2ll l34j . The differ- 
ential transmission due to the time-varying strain field 
can be written as 



/ AT 



(*) 



LA 



dz s{z, t)F(z; u>) 



(26) 



F{z\lo) is called the sensitivity junction. A schematic 
sensitivity function is shown in Fig|3| Here for simplicity, 
we assume a binary form of the sensitivity function. As 
the photoexcited acoustic waves propagate outward, its 
strain field overlaps with the sensitivity function periodi- 
cally with a time constant equal to the oscillation period. 
This periodic overlapping results in the observed oscilla- 
tions of the probe differential transmissions. Also shown 
in Fig|3 (the lower part) is the calculated probe trans- 
mission change using the sensitivity function and strain 
fields shown on left of the figure. The result agrees well 
with the experimentally measured traces. 



VI. CRYSTAL ORIENTATION EFFECT 

Here we apply the macroscopic approach developed in 
Section 4 to investigate coherent phonon generation in 
nitride-based MQWs with arbitrary growth directions. 
The current crystal-growth technique has made it pos- 
sible to grow wurtzite nitrides along orientations other 
than the conventional c-axis. Optical and electrical ex- 
perimental studies also reveal many spec ial prop erties of 
GaN with different crystal orientations [33, |3(| . Theo- 
retical studies of crystal orientation effects on wurtzite 
semiconductor band structures have been reported by 
many authors, e.g. Refill We will show in the following 
that, for MQWs with certain orientation angles, coher- 
ent transverse acoustic (TA) phonons will be excited and 
dominate the coherent LA phonon signals. This THz 
shear acoustic wave might have special applications for 
picosecond ultrasonics. 

As discussed previously, the photogenerated space 
charge density is important for the generation of coherent 
phonons. Here we first investigate the orientation effects 
on the space charges. The crystal orientation effects on 
the subband structure of strained wurtzite multiple quan- 
tum wells can be found in Refl3ll In the following calcu- 
lation, the In composition is x = 0.1, the well width L w 
is 22A, and the barrier width Lb is 90 A. As discussed 
in the preceding section, the piezoelectric coupling domi- 
nates the driving force. The piezoelectric coupling is due 
to the space-charge field created by separation of pho- 
togenerated electrons and holes, i.e. it depends on the 
magnitude of the space charge density p sc . In FigQjwe 



show the magnitude of the corresponding fundamental 
Fourier component \p S c(qo)\ versus the orientation angle 
9. The photoexcited 2D carrier density is fixed at 2 x 10 10 
cnr 2 . The magnitude of p sc (qo) has minima at 9 ~ 45° 
and 9 ~ 90°. This is because the built-in piezoelectric 
field vanishes at these angles. The maximum magnitude 
of the space charge density is at 9 = 0° , which is the con- 
ventional [0001] orientation, and a second local maximum 
occurs at 9 ~ 68°. 

After photogeneration the coherent acoustic waves 
leave the MQW region within a finite time interval, and 
a static carrier-induced strain field remains inside the 
MQW. In the following, we shall mainly consider the ra- 
diating component. We define the strain field associated 
with mode A as 



(27) 



Since the driving forces on the right hand side of Eq. Ij25(l 
satisfy the sum rule /_ f\(£,t)d£ = 0, the strain field 
s\ has a well-defined Fourier transform, which obeys the 
following driven harmonic oscillator equation 



dt 2 



n 7 



+ ut(q)s x (q,t) = ^fx(q,t), (28) 
Po 



where s\(q, t) and f\(q, t) are the Fourier transform func- 
tions of the strain field s\(£,t) and force f\(£,t). The 
oscillating frequency to\(q) = c\q. Due to the periodic 
carrier distribution within the MQW, the Fourier trans- 
form of the driving force f\(q) has peaks at q — mqo, 
where m is an integer and qo = 2ir/(L w + Lb). The 
finite linewidth of each peak f\(mqo) is due to the fi- 
nite number of wells and the inhomogeneous distribution 
of photogenerated carriers. This finite linewidth corre- 
sponds to the observed decay time constant of oscillating 
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FIG. 4: Crystal orientation dependence of the photogener- 
ated space charge density within the piezoelectric MQW. 
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traces. The solution to the above oscillator equation can 
be expressed as a superposition of displacive cosinusoidal 
oscillations 



s x (q,t) = ' / dr IX ™ j {l-cos[^(g)(t-r)]}. 

(29) 

The time derivative dfx/dr is proportional to the carrier 
generation rate, which is in turn approximately propor- 
tional to the pulse shape. Let I(t) be the normalized 
pulse shape function, i.e. I(t)dt — 1, and I(w) be 

its Fourier transform. The solution (|29[l can be approxi- 
mated as 



§ x (q,t) = S° x {q){l - I[uJ X {q)} COs[w A (?) t}}, 



(30) 




15 30 45 60 75 
Propagation Angle 6 (Degree) 



90 



(31) 



Where f\{q) = f\{q,t 3> r p ) and t p is the pulse width. 
There are two terms on the right-hand side of Ea. (l30ll . 
The first term, s°, corresponds to the static carrier- 
induced strain field [30|, This steady state field can be 
obtained by solving the following equation 



Poc\ 



(32) 



The second term on the right-hand side of Eq. jjHfll cor- 
responds to the radiating part of the strain field, i.e. the 
observed coherent phonon oscillations. The magnitude of 
the coherent acoustic wave is proportional to the steady- 
state component s ° (q) , and is reduced by the pulse factor 
I[u>\(q)]. In the following, we shall investigate the orien- 
tation effect on the magnitude of the first-order harmonic 
component, i.e. s\(qo). 

As discussed in Section 4, only two acoustic modes 
will be excited by the photogenerated carriers. The two 
modes are quasi-LA and quasi- TA modes. We shall sim- 
ply refer to them as LA and TA modes. The total acous- 
tic field can thus be expressed as 



u(£, t) = mla(£, t)w LA + u T a(£> t)wTA- 



(33) 



Where ula and -uta are the mode amplitudes of the 
LA and TA phonons, respectively. We shall first com- 
pare their modal properties versus propagation direc- 
tions. The fundamental oscillation frequencies of the 
two modes are wla = cla^o and wta = cta^Oi respec- 
tively. In FigO we show the crystal-orientation depen- 
dence of the normalized oscillating frequencies for both 
the LA (solid line) and TA (dotted line) modes. The 
normalization was made relative to the oscillation fre- 
quency of an LA mode along the [0001] direction with 
a fixed wavevector q$. The oscillation frequency of the 
TA mode is approximately half that of the LA mode. 
Both frequencies are symmetric with respect to 9 = 45°, 
while the frequency of the LA (TA) mode has a mini- 
mum (maximum) at this orientation angle. Since both 



FIG. 5: Fundamental oscillation frequencies of the LA (solid 
line) and TA (dotted line) phonons as a function of the crystal 
orientation 9 for wurtzite GaN. 



oscillating modes will contribute to the absorption modu- 
laton through piezoelectric field or deformation potential 
coupling, a beating phenomenon will be observed in the 
time-resolved pump-probe measurements. However, the 
amplitude of the beating depends also on the relative 
magnitudes of the LA and TA modes, as we shall discuss 
later. 

We now introduce the effective piezoelectric and defor- 
mation potential coupling coefficients for the two eigen- 
modes. The effective coupling coefficients appear in the 
driving force equations (A=LA, or TA) as 



fx. 



piczo 



_ Psci 
S 



(A.dcf 



E 



dp v 



(34) 



(35) 



FigEjA) shows the orientation dependence of the nor- 
malized effective piezoelectric coefficients e~LA and exA 
for the coherent LA (solid line) and TA (dotted line) 
modes, respectively. The normalization was made with 
respect to e-^x{9 = 0) = e 33 . The effective LA mode cou- 
pling constant, eLA, has a maximum at 6 = 0, i.e. the 
[0001] direction, and nodes at 6 ~ 40° and 9 = 90°. Thus 
the magnitude of the coherent LA phonon oscillations is 
expected to be small at these node angles. As for the 
TA mode, &ta has nodes at 9 = and 9 ~ 65°. Since 
the deformation potential driving force for TA modes is 
very small for all orientation angles, the TA mode can- 
not be excited in the conventional [0001] direction due 
to exA = 0. However, it is expected that in the [1010] 
direction, i.e. 9 — 90°, the generation of the TA mode 
will dominate over that of the LA mode. In FigH^B), we 
show the normalized effective coefficients for electron de- 
formation coupling, d e la and (i e .TA, as functions of the 
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FIG. 6: Normalized effective piezoelectric (A) and deforma- 
tion potential (B) coupling coefficients of LA (solid line) and 
TA (dotted line) phonons as a function of the crystal orien- 
tation 9 for wurtzite MQW. 



propagation angles. The normalization was again made 
with respect to the d e ^A of the LA mode in the [0001] 
direction. As discussed previously, the deformation coef- 
ficient is approximately proportional to the projection of 
the propagation unit vector n on the mode polarization 
vector w. It can be seen from Fig||jB) that the normal- 
ized coefficient of the LA mode is almost equal to 1 for all 
orientations, while that of the TA modes is almost zero. 

Now we shall investigate the orientation dependence of 
the driving forces for the LA and TA modes. We consider 
a MQW with 10% In composition, a well width L w = 22 
A, and a barrier width Lb = 90 A. We assume a constant 
2D photogenerated carrier density of 2 x 10 10 cm -2 for all 
orientations. We can thus compare the intrinsic effects of 
crystal orientation on the generation of acoustic phonons. 
Please note that the optical absorption properties depend 
on the orientation, one needs different pump intensities 
to achieve the same carrier density at different orienta- 
tion angles. The orientation dependence of the photo- 
generated carrier densities and space charge density cor- 
responding to the same MQW parameters are discussed 
above, e.g. Fig0]for the space charge density. Those re- 
sults combined with the effective coefficients are now used 



to derive the driving forces for each mode. In FigjJfA), 
we first show the magnitude of the driving forces ver- 
sus orientation angle for the LA mode. Here we only 
show the Fourier component corresponding to the funda- 
mental wavevector go- As can be seen from the figure, 
the piezoelectric force /LA,piezo dominates the generation 
of coherent phonons and the deformation potential cou- 
pling force /LA.dof is almost independent on the orienta- 
tion angles. From Ea. (|34(l . the piezoelectric driving force 
/LA.piczo is proportional to the effective piezoelectric co- 
efficient and the space charge density. The magnitude of 
the space charge component has a minimum near 8 ~ 45° 
[See FigQ] and the effective LA mode piezoelectric co- 
efficient has nodes at both 8 ~ 45° and 8 = 90° [see 
FigElA)]. As a result, the piezoelectric force for the LA 
mode vanishes at both orientation angles. However, con- 
trary to the case of the TA modes (discussed below), the 
LA mode has a constant contribution from deformation 
potential coupling for all orientation angles. 

FigEjB) shows the magnitudes of the driving forces 
versus orientation angles for the TA mode. The defor- 
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FIG. 7: Driving forces as a function of the crystal orientation 
9 for (A) LA and (B) TA phonons. 
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FIG. 8: Magnitudes of the generated coherent LA (cirle) and 
TA (triangle) phonons as a function of the crystal orientation 
9 for a fixed photogenerated 2D carrier density. 



mation coupling force /xA.def is vei T small compared to 
the piezoelectric force /TA.piozo and the piezoelectric force 
vanishes at angles corresponding to the nodes of the ef- 
fective TA piezoelectric coefficient. A kink-like minimum 
at 9 ~ 45° in the magnitude of the piezoelectric force is 
due to a minimum in the space charge component. The 
piezoelectric force is largest at angles ~ 30°. 

Finally, in Fig|H| we show the resultant strain ampli- 
tudes of both modes versus the orientation angles. It 
can be seen that although only the LA mode is present 
in the [0001] direction, for other orientations both LA 
and TA modes will be excited. The generation of the 
TA mode even dominates over that of the LA mode for 
certain orientation angles, especially at 9 ~ 30° and 90° 
(the [1010] direction). In addition to the well-known LA 
acoustic wave generated along [0001] direction, our re- 
sults suggest that shear acoustic waves can be excited 
in wurtzite semiconductor MQW's with suitably chosen 
crystal orientations. For example, TA-dominant acoustic 
waves can be generated by photoexciting a (1010)-grown 
MQW. These findings may have potential applications in 
picosecond ultrasonics. 



femtosecond pulse excitation, the large built-in piezoelec- 
tric field in the InGaN well region spatially separates 
the photogenerated electrons and holes. This separation 
causes a space charge density that partially screens the 
original electric field in the MQW. Due to piezoelectric 
coupling, the corresponding equilibrium lattice state is 
shifted and a displacive coherent LA phonon oscillation 
is then initiated. The photogenerated carrier distribu- 
tion is confined within the well region and is a periodic 
function with wave vector go = 2tt/(L w + Lb), where 
L w and Lb are the well and barrier widths, respectively. 
Coherent phonon oscillation, with the selected acoustic 
phonon mode corresponding to the wave vector qo, can 
then be initiated. The phonon has an oscillation fre- 
quency luq = c s qo. The induced acoustic phonon oscil- 
lation results in piezoelectric field modulation and then 
causes absorption variation through the quantum con- 
fined Franz-Keldysh effect. The observed decay time con- 
stant of the measured trace indicates the time required 
for the generated coherent phonons to leave the MQW 
region. 

We also reviewed the recently proposed microscopic 
theory for the generation and propagation of coherent 
LA phonons in wurtzite semiconductor MQW's. The co- 
herent acoustic phonon oscillation is a macroscopic co- 
herent population of the phonon states. In this respect, 
it is similar to a laser photon field, which is a macro- 
scopic coherent state of photons. Motivated by appli- 
cations to more complicated coherent acoustic devices, 
we developed a macroscopic theory for the generation 
and dynamics of coherent acoustic phonons in wurtzite 
MQW's. The approach is based on macroscopic con- 
tinuum constitution equations for piezoelectric wurtzite 
semiconductors. Starting from Poisson equation and the 
dynamic elastic equation, a vector loaded wave equation 
was obtained. By projecting the corresponding equa- 
tion to eigenvectors of the elastic Christoffel equation, 
the loaded string equation is derived. Only two acoustic 
eigenmodes can be excited by photogeneration. They are 
the quasi-longitudinal and quasi-transverse modes. We 
also applied the macroscopic approach to study crystal- 
orientation effects on the generation of coherent acoustic 
phonons. It is found that although only the LA phonon 
is photoexcited as the crystal is grown along the [0001] 
direction, at other orientations, e.g. [1010], the genera- 
tion of TA phonon is favored. This coherent transverse 
wave may have potential applications in THz ultrasonics. 



VII. CONCLUSION 

In this paper we have provided an introductory review 
of the experimental and theoretical aspects of coherent 
acoustic phonon generation in nitride-based semiconduc- 
tor heterostructures, with particular application to In- 
GaN/GaN MQWs. 

We first presented the observation of coherent 
LA phonon oscillations in InGaN/GaN MQWs from 
transmission-type pump-probe measurements. With UV 
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